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The dyadic Green’s function of the inhomogeneous vector Helmholtz equation describes the field 
pattern of a single frequency point source. It appears in the mathematical description of many areas 
of electromagnetism and optics including both classical and quantum, linear and nonlinear optics, 
dispersion forces (such as the Casimir and Casimir-Polder forces) and in the dynamics of trapped 
atoms and molecules. Here, we compute the Green’s function for a layered topological insulator. 
Via the magnetoelectric effect, topological insulators are able to mix the electric, E, and magnetic 
induction, B, fields and, hence, one finds that the TE and TM polarizations mix on reflection 
from/transmission through an interface. This leads to novel field patterns close to the surface of a 
topological insulator. 

PACS numbers: 78.20.-e, 78.20.Ek, 78.67.Pt, 42.25.Gy 


I. INTRODUCTION 

Topological insulators are a class of time-reversal sym¬ 
metric materials that display non-trivial topological or¬ 
der and are characterized by an insulating bulk with pro¬ 
tected conducting edge states mm- This type of mate¬ 
rial was first predicted [3] and then observed |1] in 2D 
in HgTe/CdTe quantum wells and subsequently in 5D 
in Group V and Group V/VI alloys that display strong 
enough spin orbit coupling to induce band inversion - 
Bii_xSbx in the first instance mm and then in Bi 2 Se 3 , 
Bi 2 Te 3 and Sb 2 Te 3 [3 [S] to name but a few examples. 
Owing to their unusual band structure, these materi¬ 
als display a number of unique electronic properties, the 
most notable of which is the quantum spin hall effect 
where quantized surface spin currents are observed even 
though the usual charge currents are absent PIIQ]. 

In addition to their interesting electronic properties, 
topological insulators also display a number of unusual 
electromagnetic properties. Specifically, topological in¬ 
sulators have the ability to mix electric, E, and mag¬ 
netic induction, B, fields [Hill], a feature which has a 
pronounced affect on the optical response of the mate¬ 
rial (TSHlS]. In particular, this magnetoelectric E — B 
mixing allows one to realise an axionic material [Mini- 
Such materials are described by the Lagrangian density 
>Co + Caxion, where Co is the usual electromagnetic La¬ 
grangian density and Caxion is a term that couples the 
electric and magnetic induction fields. This additional 
electromagnetic interaction reads 

Caxion = -^E(r,a;) • B(r,w), (I) 

[lOC 
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where a is the fine structure constant and 0(r,a;) is 
termed the axion field in particle physics (although, as 
far as electromagnetism is concerned, it merely acts as 
a space and frequency dependent coupling parameter). 
In order to realise such a material in a topological in¬ 
sulator, a time symmetry breaking perturbation of suffi¬ 
cient size must be introduced to the surface to induce a 
gap, thereby converting the material into a full insulator. 
Such a gap can be opened by introducing ferromagnetic 
dopants to the surface (12% Fe doping in Bi 2 Se 3 leads 
to a mid-infrared gap of 50meV/25p,m [13]) or by the 
application of an external static magnetic field [19] ■ In 
such a time-reversal-symmetry-broken topological insula¬ 
tor (TSB-TI) the constitutive relations are altered and, 
hence, the optical properties of the material change dra¬ 
matically. 

Here, we derive the electromagnetic Green’s function 
for a layered TSB-TI. The electromagnetic Green’s func¬ 
tion is the solution to the vector Helmholtz equation for 
a single frequency point source and can be used to gen¬ 
erate general field solutions for an arbitrary distribution 
of sources. This function has a wide range of applica¬ 
tions in both classical [MIH] and quantum optics [33U23] 
and is an important component in studies of linear [26j 
and non-linear [13 HH] optics, Gasimir [39] and Casmir- 
Polder [MEI] forces, decoherence [H] and the dynam¬ 
ics of trapped atoms [33] and molecules [MjES]- Thus, 
knowledge of the Green’s function is of value to a great 
many fields. 


II. MAXWELL EQUATIONS 

As with all electromagnetic studies, we begin with the 
Maxwell equations and constitutive relations for the ma- 
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terial in question. For a TSB-TI these are [niiiiin] 


V-B(r,w)=0, (2) 

V X E(r, w) — za;B(r, w) = 0, (3) 

V • D(r,a;) = p(r,w), (4) 

V X H(r,w) + ia;D(r, w) = J(r,a;), (5) 


D(r,a;) = £o£(i’, w)E(r, w) 
a 0(r, (jj) 


H(r,o;) = 


TT fJ,oC 

B(r,w) 


B(r,a;)+PAr(r,a;), (6) 


a 0(r, oj) 

TT floC 


E(r,w) - MAr(r,a;), (7) 


where a is the fine structure constant and e(r, w), /r(r, uj) 
and 0(r,a;) are the dielectric permittivity, magnetic 
permeability and axion coupling respectively, the lat¬ 
ter of which takes even multiples of tt in a conventional 
magneto-dielectric and odd multiples of tt in TSB-TI, 
with the magnitude and sign of the multiple given by 
the strength and direction of the time symmetry break¬ 
ing perturbation. The PAr(r, w) and MAr(r, w) terms are 
the noise polarization and magnetization, respectively. 
These terms are Langevin noise terms that model ab¬ 
sorption within the material [53] . These relations can be 
derived from the Lagrangian density in Eq. Q [T3]. Us¬ 
ing the above constitutive relations, one can show that 
the frequency components of the electric field obey the 
inhomogeneous Helmholtz equation 


V X —- -V X E(r,a;) - ^^(r, a;)E(r, w) 

— i—— [V0(r,a;) x E(r,a;)] 

C TT 

= iw/xo [J£;(r,a;)-1-Jjv(r,w)], (8) 


where J£;(r,a;) is the source term for electromagnetic 
waves generated by external currents and J^vlr, w) = 
—iojPN{r,uj) -I- V X Mjv(r,a;) is the source term for 
electromagnetic waves generated by noise fluctuations 
within the material. If the axion coupling is homoge¬ 
neous, 0(r,a;) = Q{oj), then the last term on the left- 
hand side vanishes and one finds that the propagation of 
the electric field is the same as in a conventional magneto¬ 
dielectric. As a result, electromagnetic waves propagat¬ 
ing within a homogeneous TSB-TI retain there usual 
properties - dispersion is linear, the phase and group ve¬ 
locities are proportional to the usual refractive index, the 
fields are transverse and orthogonal polarizations do not 
mix. Thus, the effects of the axion coupling are only felt 
when the axion coupling varies in space. For layered, ho¬ 
mogeneous media this will occur only at the interfaces 
where the properties of the medium change. 


III. FRESNEL COEFFICIENTS 

An important set of functions for any layered media 
are the Fresnel coefficients for reflection and transmis¬ 
sion at each interface. These functions are required to 
construct the Green’s function. In fact, the Fresnel co¬ 
efficients for TSB-TIs have been studied before [HUT], 
however, the standard expression for the Green’s function 
requires a slightly different form for the coefficients com¬ 
pared to those in previous work [501IM]- Furthermore, 
certain aspects of TSB-TIs mean that the usual method 
of computing this form of the coefficients leads to incor¬ 
rect results. For these reasons, it is worth revisiting the 
derivation in some detail. 

In the derivation of the Fresnel coefficients for a con¬ 
ventional magneto-dielectric one usually defines two po¬ 
larizations; the TE polarization, where the electric field, 
E, is parallel to the interface, and the TM polarization, 
where the magnetic field, H, is parallel to the interface. 
Since, from Eq. the electric field propagation is unaf¬ 
fected by a homogeneous axion coupling one can see that 
the TE polarization is unchanged. However, from Eq. 
0, one can see that the magnetic field, H, is no longer 
perpendicular to the electric field, E. Thus defining the 
TM polarization in terms of the magnetic field, H, leads 
to two polarizations that are not orthogonal and hence 
incorrect expressions for the Fresnel coefficients. Fur¬ 
thermore, we would expect the two polarizations to mix 
at the interface via the magnetoelectric coupling, hence 
defining the two polarizations in terms of different fields 
leads to awkward expressions. The simplest approach is 
to work solely with the electric field, E. Thus, for media 
layered in the z direction and light incident in the x — z 
plane, we define the TE polarization as the polarization 
with Ey ^ 0 and E^ = 0, E^ = 0 and the TM polariza¬ 
tion as the polarization with Ey = 0 and E^ 0, E^ 0, 
[See Fig. [^. 

We proceed by considering waves of a specific {TE, 
TM) polarization incident on an interface between two 
homogeneous isotropic TSB-TIs. By matching the waves 



FIG. 1: The interface between two topological insulators. 
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in each half-space using the electromagnetic jump condi¬ 
tions 


z X El = z X E 2 , (9) 

z X Hi = z X H 2 , (10) 


which relate the transverse components of the electric 
and magnetic fields on either side of the interface, the 
Fresnel coefficients can be found. 

First we consider a TE polarized plane wave incident 
on the interface from layer 1 [See Fig. [^. The electric 
field ansatz for each region is 


tions for the fields at the interface at z = 0 


1 -I- Rte,te = Tte,te, 


- ^TM,TE — 

ni 712 


2,2 r 


f-TM,TE^ 


—— [1 — Rte.te] H— Oi^^Rtm,te 

fll TT Til 

kz^2rp 

= - J^TE,TE -*32- J-TM,TE, 

^2 E 712 

— Rtm.te -©1 [1 + Rte,te] 

Ml 7^ 

= —Ttm,TE - Q2TtE,TE, 

M2 77 


(21) 

( 22 ) 


(23) 


(24) 


E.,i = (11) 

fci 

Ey^i = -Eo , (12) 

E.,1 = (13) 

ki 

E.,2 = (14) 

K2 

Ey,2 = (15) 

E.,2 = Eo’^e-^’^^-^^+^'^-^TTM,TE, (16) 

k2 


where the k^/ki = cos 4>r, kp/ki = sin^^, fcz /^2 = cosc^t, 
kp/k 2 = sin (j)f It has been previously shown that Snell’s 
law holds for TSB-TFs so (j)i = (j)r [HI HI]. From Eqs. 
^ and Q we obtain 


I = -Eo^^ l^-^k.,^z+ik^x _ ^^1 

/ioMiw 


Ct ©1 fczp ik^,iz+ik„x p 

+ i^o- ^ —e ’ p Remote, 

TT //qc ki 

Hy,i = -Eo— 

/ioMiw 

©1 ^^—ikz,iz+ikp 

\ flQC 


(17) 


+ Eo -|e ■--p- +e 


X _|_ ^ikz,iz-\-ikpX 


JET _ ^Z,2 —ikz,lZ-\-ikpXrp 

Rx,2——Ro -e ’ P J-TE,TE 


/i0^2W 

Q: ©2 kz^2 —ikz iz+ik 


— Eq - j^e 

TT /IqC k2 

h 




TM,TE, 


TE,TE\ ) 

(18) 


(19) 


JLf _ TP '2 —ikzAZ + ikpXrp 

Ey,2 — -Eq -e • P J^TM,TE 


M 0 M 2 W 

©^ —ikz.iz+ikpQ 


+ Eq -e 

TT ^oC 


T, 


TE,TE- 


( 20 ) 


where we have used the dispersion relation k = nLojc. 
Solving the above system of equations gives 


Rte,te = 
Remote = 
Tte,te = 
Ttm,te = 




2M2^z,llle 


2^2712^2,1^ 


(M2fcz.l + k-lkz,2)^e + fcz,lfcz.2A2 


:z.2A2’ 

(25) 

:z.2A2’ 

(26) 

:z.2A2’ 

(27) 

:z.2A2’ 

(28) 

MiM2(fcz,ie2 + 


kz, 2 £i): with the factors of e appearing via the defi¬ 
nition of the refractive index, = jis. Note that 
when ©2 — ©1 —>■ 0 (i.e when the axion couplings van¬ 
ish or when they are the same across the interface), 
Rtm,tE:Txm,te 0 and Rte,te and Txe,te reduce 
to the usual reflection coefficients for normal magneto¬ 
electric materials |20l |24| . 


Next, we consider a TM polarized plane wave incident 
on the interface from layer 1 [See Fig. [^. The electric 
field ansatz for each region is now 


Ex,i 

Ey,l 

Ez,i 

Ex,2 
Ey,2 
Ez,2 


k 1 

_ [ —ikz,iz-\-ikpX AkzTZ+ikpX td 1 

— Eq-^ [e • p — e • p Etm,tm\ , 

Ki 



(29) 

_E^Ek^A^+^kpXR^E^TM, 

(30) 

ki 

1 

(31) 

TP ik 2Z-\-ikpXrp 

Rq—j -e ’ ^ J-TM,TM, 

k2 

(32) 


(33) 

k 

77' P ^ — ikzOZ+ikpXrp 

k2 

(34) 


From Eqs. ([^ and (101 we can find the boundary condi- 


where, again, components of the wavenumber are related 
to the angles of incidence, reflection and transmission and 
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Hy^l — —Ecl 


ki 


fioHiUJ 

a 01 


y—ikz iz-\-ikpX I ^ikz iz-\-ikpX 


Snells law holds. From Eqs. ^ and Eq. Q we obtain 

CX 01 fcz^l r ij. 2 +zfepa; „ik^AZ+ik„x-n 1 

— iio-;— e ' p — e ’ ^ Etmtm , 

TT /ipC fci 

(35) 

(36) 

(37) 

(38) 


+ Eo- — 


TT /ioC 


TE,TM, 


ZJ _ ZP ^ — 'ikz.2^-\-‘i‘kr,Xrj^ 

tlx,2 — —^0 . € J~TE,TM 


a 02 22+*fcp 


- Eo- 

TT /XqC fc2 




TM,TM, 


ZJ _ Z7’ ^ ^—'i-kzAX+iknXji 

My 2 — —Mo -e • P xtm,TM 

M 0 M 2 W 

Q! 02 i}; z+ikj,xrp 

+ Mo -e • p 1te,tm- 

TT fJoc 


From Eqs. § and ( | 10 [ ) we can, again, find the boundary 
conditions for the fields at the interface at ^ = 0 


[1 — Rtm,tm] = -^Ttm,tm, 

ni n2 

Rte,tm = Tte,tm, 

p “o ri p 1 

- MtE,TM -tJi- i — Utm,TM\ 

III TT ni 

kz,2rp Of fcz,2rp 

—- Xte,TM -C'2- xtm,TM, 

fJ.2 TT n2 

— [1 + T^tm.tm]- QiRte,tm 

Ml TT 

= TtmtM -027tE,TM, 

M2 TT 


(39) 

(40) 

(41) 

(42) 


where, once more, k = niojc has been used. Solving the 
above system of equations gives 


the transmission coefficient differs from the above result 
by the ratio of the impedances of the two layers. This is 
because the TM coefficients are derived using the H-field 
instead of the E-field.) 

The Fresnel coefficients for the energy flux can be 
found by comparing the z-component of the Poynting 
vector, S = E X H, on each side of the interface. One 
finds that they are related to the above field coefficients 
via 




j I > 

, _ fcz,2 Ml Irp |2 

^i,j — 1 ,, 1 

Kz,l M2 


(47) 

(48) 


where i,j € TE,TM and the prefactor in the transmis¬ 
sion coefficients accounting for the change in flux area 
as the field passes through the interface. Figure shows 
the reflection, and transmission, as a function of 
incident angle for a 600nm optical plane wave, incident 
from the vacuum, encountering a conventional magneto¬ 
dielectric with £ = 16 and M = 1 (such values are simi¬ 
lar to those for 612803 at high frequencies [SB])- In this 



Rtm,tm 

Rte,tm 

Ttm,tm 

Tte,tm 


{£2kz,i 

- £lkz,2)^y + k:z,lkz,2^'^ 

i£2kz,i 

+ £lkz,2)^iz + kz,lkz,2PP^ ’ 


-2n2nikz,ikz,2P^ 

i£2kz.i 

+ sikz, 2)^11 + fcz,iA:z,2A2 ’ 

n2 

‘^£lkz,l^y 

ni (£2A 

■2,1 + £lkz, 2)^11 + fcz,lfc2,2A^ 


-2n2nikz,ikz,2Pk 


{£2kz,l + £lkz,2)^y + kz,lkz,2PP ’ 


(43) 

(44) 

(45) 

(46) 


where A = aMiM 2(02 — 0i)/7’‘ and 11^ = MiM 2 (fcz,iM 2 + 
kz, 2 k'i), with the factors of £, again, appearing via the 
definition of the refractive index. Once again, when 
02 — 01 —>■ 0 (i.e when the axion couplings van¬ 
ish or when they are the same across the interface), 
Rte,tmiTte,tm —t 0 and Rtm,tm and Ttm,tm reduce 
to the usual reflection coefficients for normal magneto¬ 
electric materials [IHl 121] ■ (Note that in [501 1 unlike |53| , 


FIG. 2: (Color online) The % reflection (solid) and transmis¬ 
sion (dashed) for TE (red) and TM (blue) polarized waves at 
a vacuum- (layer 1) magneto-dielectic- (layer 2) interface as a 
function of the incident angle, 4>i- Here, pi = p 2 = 1, ei = 1 
and £2 = 16. 


case the mixing coefficients vanish and the polarization 
state of the incident light is preserved by the interface. 
It is easy to see that for incident TE polarized light 
fTE,TE + tTE,TE = 1 hence the TE energy flux is pre¬ 
served at the interface (a similar expression holds for T M 
polarized light). In comparison. Fig. shows the reflec¬ 
tion, TiA, transmission, and mixing, rij/Tj (i j), 
as a function of incident angle for a plane wave of simi¬ 
lar wavelength, incident from the vacuum, encountering 
a TSB-TI with £ = 16, p = 1 uud 02 = tt. In this 
case the mixing coefficients are non-zero and the polar¬ 
ization states of the incident light mix at the interface. 
Finally, one can show that for incident TE polarized 
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FIG. 3: (Color online) (a) The % reflection (solid) and trans¬ 
mission (dashed) for TE (red) and TM (blue) polarized waves 
at a vacuum- (layer 1) TRSB-TI- (layer 2) interface as a func¬ 
tion of the incident angle, (pi. (b) The % reflection (solid) 
and transmission (dashed) for TE ^ TM mixing (red) and 
TM ^ TE mixing (blue) at a vacuum- (layer 1) TRSB-TI- 
(layer 2) interface as a function of the incident angle, (pi. In 
both cases, /ri = p2 = 1, £i = 1, £2 = 16, 0 i = 0 and 02 = tv . 


light, tte.te + fTM.TE + tTE,TE + txM.TE = 1 hence 
TE energy flux is still preserved at the interface (again 
a similar expression holds for TM polarized light). 

In order to better understand the mixing coefficients, 
it is informative to look at the case of a pure TSB-TI 
where the permittivity and permeability are that of the 
vacuum and only the axion coupling changes on the in¬ 
terface. This allows one to remove the magneto-dielectric 
effects from the system and isolate the effect of the axion 
coupling. In this case kz,i = kz ^2 and the reflection and 
transmission coefficients reduce to 

-A2 

Rte,TE = —Rtm.TM = y——yW) (49) 

4 -b 

4 

Tte,te = Ttm.tm = ■ , A 9 ^ (^ 0 ) 

4 -b A^ 


RtM,TE = RtE,TM = —Ttm,TE = Tte,TM = . , . n ■ 

4 -b 

(51) 

(Similar expressions were found in |17].) One can see 
that in the pure TSB-TI limit the reflection and trans¬ 
mission coefficients are no longer a function of incident 
angle and, hence, the angular dependence of the coef¬ 
ficients is a result of the magneto-dielectric properties 
of the material rather than the axionic properties. As 
A « a, the energy flux reflection and transmission coeffi¬ 
cients are « a‘^/16 « 10“^° and « 1 respectively. 
Thus one sees near perfect transmission. However, since 
the axion coupling changes on the interface one still has 
mixing, the magnitude of which is equal in transmission 
and reflection Vij = tij « 0^/4 « 10“^. 

As a slight diversion we briefly consider reflection, 
transmission and mixing for large values of A. Although 
such large values are probably not realizable with a topo¬ 
logical insulator, this limit is useful in understanding the 
affect of the axion coupling and may bear some relation 
to treatments of the axion coupling as a fundamental 
field. Figure 1^ shows the reflective, transmissive, 
ti^i, and mixing, r^j-, /tij, coefficients for the expressions 
in Eqs. ( [4^ - ( [sTj ) as a function of A. For vanishing A 
one sees that the reflection and mixing coefficients van¬ 
ish and one has perfect transmission. For A —> 00, the 
transmission and mixing coefficients vanish and one ap¬ 
proaches a perfect mirror. Thus, for large changes in the 
axion coupling, the interface becomes purely reflective 
and no mixing occurs. Further, one sees that the maxi¬ 
mum mixing occurs when the reflection and transmission 
coefficients are equal. However, if one changes the rela¬ 
tive impedances of the layers, this maximum is shifted. 
Increasing the impedance (increasing /r relative to e) in¬ 
creases the mixing and shifts the peak to A values lower 
than the reflection-transmission crossing point [See Fig. 
(a)], while lowering the impedance (increasing e rela- 



FIG. 4: (Color online) The % reflection (red) and transmis¬ 
sion (blue) and mixing (purple) as a function of A for a pure 
TI, with Pi = P 2 = £1 = £2 = 1. 
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FIG. 5: (Color online) The % reflection (red), transmis¬ 
sion (blue), reflective mixing (black) and transmissive mix¬ 
ing (purple) as a function of A for a TSB-TI, with (a) 
= £i = £2 = 1 and fi 2 — 2 and (b) /ri = /r 2 = £2 = 1 
and £2 = 2 for normal incidence angle (hence TE and TM 
polarizations are indistinguishable). 


IV. THE GENERALIZED FRESNEL 
COEFFICIENTS FOR MULTILAYERED MEDIA 


To find the reflection and transmission coefficients for 
multi-layered media we follow the method of Ref. [20] 
and first consider a three layered medium. The wave in 
layer 1 can be written, in the TE/TM basis as 


El = 


g-ifej.izj _|_ gife2,i(z-2di2)^ 


■E, 


0 . 1 ) 


(54) 


where / is the unit matrix and R ^2 is the generalized 
reflection matrix at the layer 1 - layer 2 interface, whose 
form is unknown. Similarly in layer 2 we have 


Eo = 


g-ifez,2ZJ _j_ gjfe2,l(z-2d23) 


23 


■E, 


0 . 2 ) 


(55) 


and in layer 3 we have 

E3 = A3 • • Eo, 3. (56) 

By considering upward and downward propagating waves 
in each layer one can compute the unknown matrices. 
The downward propagating wave in layer 2 is a conse¬ 
quence of the transmitted wave from layer 1 and the re¬ 
flected wave from layer 2 - layer 3 interface 


g-i/c..2<ii2E^ 2 = ■ Eo.i 

-b • R 23 ■ Eo.2) (57) 

which can be solved for Eo,2 to give 

Eo ,2 = ■ T ,2 ■ Eo.i, (58) 

where 

M 2123 = I - • ^ 23 ) ( 59 ) 

and the power of —1 implies the matrix inverse. The 
upward propagating wave in layer 1 is a combination of 


five to fj.) leads to a decrease in the mixing and shifts the 
peak to A values larger than the reflection-transmission 
crossing point [See Fig. (b)]. Thus, we see the mixing 
is enhanced by the magnetic response of the material and 
suppressed by the electric response. 

Finally, it will be convenient for the rest of this study 
to write the Fresnel coefficients in matrix form 


E = 

T = 


Rte,te 

Rtm,te 


Rte,tm 

Rtm,tm 


Tte,te Tte,tm 
Ttm,te Ttm,tm 


(52) 

(53) 


These matrix transformations act on the field vector E 
whose components refer to the TE and TM polarizations 
respectively. 
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reflected waves from the layer 1 - layer 2 interface and 
transmitted waves from layer 2 


— ifc^ i(ii2 


Sl12 


Eo.i = • Eo,i 

I ^ikz,2{di2 — ‘^d22.)'T' T? 1 ? 

+ e J^21 ■ ^23 -■^0,2, 


(60) 


Adding further layers below layer 3 merely requires one 
to replace ^23 with ^23- Thus, one obtains a recursive 
relation for the reflection coefficient 


which with the help of Eq. (581 can be used to solve for 
the generalized reflection coefficient 


if — P _|_„-2ifci,2(d23-<ii2)T-’ .R -T 

i!:12 — iil2we ^21 i!:23 iii-2123 ^12' 


(61) 


R 


= R, 




_j_ g-2ifc^_i+i(cii+i,i+2-<ii,i+l) 


R, 


—2 + 1,2 —2 + 1 , 2+2 


M 


-1 

•2 + 1,2,2+1,2+2 ’^2,2+l5 


T 


(62) 


and hence one can generate the reflection coefficient for a medium with any number of layers. For transmission, from 
Eq. (58), downward going waves in layer 2 are given by 


e-*fc..2A2Eo 2 = e-*'=-i''-M2+23 ■ 1:12 • Eo.i, (63) 

and similarly for the following layer 

„ —*^2,3^23 1,' _ *fcz,2d23 A/f~l T' TP _ ifcz,2(ci23 —di2)„ —jfcz,idi2 A/t~^ T A/f~^ T TP fRA\ 

e Ji'O.S — e 2^3234 • A23 -*^0,2 — e 'e • 1^3234 • A23 ^^2123 A12 lD4j 

Thus, one can see that the generalized transmission coefficient reads 

N 


qh _ „ifcz,iv(div,N+i —div-i,jv) TT *fcz,i(di,i+l —di_l i) —1 rp 

zS-lAT — 6 


(65) 


r 


V. EMBEDDED SOURCES 




d+ 

d. 


/source (TE, TM) 


Zk 




R. T 




FIG. 7: Layered media with an embedded source. 

Consider a source embedded in a layered media. The 
source produces a unit amplitude electric held with polar¬ 
ization, Eq, in the TE/TM basis. The general expression 
for the held E in the same layer is given by 


E = F-Eo = 


ikzZ 


B 


^ikzZ : 


D -Eq. (66) 


The terms with coefficients B correspond to downward 
propagating waves and is a result of rehections from the 
surface at d+ and the terms with coefficients D corre¬ 
spond to upward propagating waves and is a result of 
rehections from the surface at d_. Thus, at the upper 
interface, we have 


^g-ifc,d+ , [jgifez|d+-^'| +^gifezd+ 


(67) 


where ,R_|_ is the generalized rehection matrix at the 
surface and at the lower interface, we have 


De 


,ikzd_ 


= R. 


J^ik^\d_-z'\ ^^-ik^d_ 


( 68 ) 


where R_ is the rehection matrix at the d_ surface. Solv¬ 
ing for B and D gives 


Be-^k.d+ ^ 


+- 


^ik^\d^-z'\^ 




R- 


_j_gifc^ -d_) gifc^ I d_-z'I ^ 

and 

^ g * fczd _ ^ 

_|_gifcz(d+-d_)gifc,|d+-z'|^ . 


, (69) 


gZfej|d_-z'|^ 


, (70) 






















with the multiple reflection coefficient reading 


—ij 




(71) 


Again, the power of —1 denotes the matrix inverse. 
Note that, unlike standard magneto-dielectrics, in gen¬ 
eral ^ _Substituting the expressions for B 

and D back into the expression for the field, noting that 
z, z' > d~ and z, z' < d+and using the definition in Eq. 
0 leads to 


F^>z'(z,z') = 


^ikzZ 


M 


-I- e 


--t 


—ikz 




'-M 


■+- 


g-ifcj(z'-2d+)^^ 


_^^ik,{z'+2d+-2d.}^^ ^ 


(72) 


for z> z' and 
F^<z/(z,z') = 


— ikz 


'-M 


+- 


-k 




-+ 


^ikz{z'-2d^) ^ 


j^^-ik,{z'+2d_-2d+)^_ 


(73) 


for z < z'. 

One can also find the electric field in a different layer 
from the source by considering the transmitted fields. 
The general expression for the field in layer n as a result 
if a source in layer m < n is 


E = F E„ = 






■E„. (74) 


The upward going field in layer n at the d„_ interface 
can be written as 


Fn — e 


— ikrriyzdn 
+ 

^ — ikm zdfi- 


'E-mn ' Fm 

)E„- •E„+-F„ 

^E.n-n+ ■ i-mn ' Fm, 


(75) 


where the field F^ at the dm+ interface is given from Eq. 
(P by 


E^ = 

„—ik^ 


■m—m+ 

jr_^gzfc„„(z'-2d„_)fo 


Fo. (76) 

Thus, the F matrix for the field in layer n is given by 


Fz>^'(z,z') = 


Ak^ 


,,zZj _|_ g-i/c„,^(2-2d„+)^ 


— ikm zdn 


M 


•n—n+ —mn 


■M 




km,zdm+ 


z-^k^,zz'J^^^k^Ad-2dz^-]fl 


■ (77) 


Similarly, the general expression for the field in layer n 
as a result if a source in layer m > n is 


E = F • E„ = 


— ik.n 


_|_ gik„,z{z-2d^^) ^ 


■E„. (78) 


The downward going field in layer n at the dn-\- interface 
can be written as 


n+'T' . 17 

=L-mn ^rn 


E„ = 

-k , • R„_ ■ E„ 


^n+ =^n- 

_ pikrr.,zd„+ ^ nn . xp 

^ ii^n+n— —mn ^m't 


(79) 


where the field E^ at the dm+ interface is given from Eq. 
(I73| by 


f' — p-ik^,zdrr,+ n> 

J-irn. — C- iy± 


•m+m— 


Akm 2-2 




—m+ 


En. 


Thus, the F matrix for the field in layer n is given by 


F^<^/(z,z') = 


— ik-n 


_|_ gifcn,z(z-2d„_)^ 


Akm zdn 


rp -ikm,zd^ + 

^ ±l±_n+n— —mn m-\-m— ^ 


Akm z 2' 


'j_^g-*fc^,.0'-2d,„+)^ 


±^m-\- 


■ (80) 


VI. THE MULTILAYERED GREEN’S 
FUNCTION 

The Green’s function is the solution to the wave equa¬ 
tion, Eq. Q, for a single frequency point source. For a 
homogeneous axionic coupling this wave equation reduces 
to 

1 A 

V X —- -V X E(r,a;)-^e(r,w)E(r,w) = zw/xoJ(r), 

/r(r,w) 

(81) 

where J(r) = j£;(r,tLi) -k Jjv(i", w) is the total current 
source with both external and noise contributions. Thus, 
the Green’s function is defined by 


V X / V X G(r,r',a;) 




-^£(r,w)G(r,r',w) =d(r-r'). (82) 

Knowledge of the Green’s function allows one to compute 
the electric field at any point from an arbitrary distribu¬ 
tion of current sources via 

E(r,a;) = J d^r'G{r,r',uj) ' (83) 


As can be seen from Eq. (811, the wave equation for 


a homogenous axionic coupling is just the usual wave 
equation for a traditional magneto-dielectric. Hence the 
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Green’s function is identical to the standard magneto- 
dielectric electric Green’s function, which, in its singu¬ 
larity extracted form, reads 


G(r,r',a;) = ^ / 


m*(r) 0 m(r') 


h h‘^ 

i^z^p 


n*(r) 0 n(r^) 


Jc ^2 
rvz^p 


gikp-(rp-rp)gifc^|2-0'| 


(84) 


For a field point at z > 0 there are two contributions, 
one from direct propagation from the source to the field 
point, which is given by the free space Green’s function 
Go(r,r' , w), and one from reflections from the surface, 
which is given by the reflective part of the Green’s func¬ 
tion i?(r,r',a;). For z < 0 the only contribution is from 
transmission at the surface, which is given by the trans¬ 
missive part of the Greens function T{r,r',uj). Thus, we 
can split the Green’s function into 3 parts 


where m(r) and n(r) are the dyadic operators 

m(r) = iVr X z, (85) 

n(r) = y Vr X Vr X z, (86) 

k 

that generate the solenoidal vector wave functions [5D], 
which are equivalent to the polarization vectors in |24j . 
Here, k is the wavevector of the wave with kp = kxX+kyy 

and kz = ^kP' — k'^. Similarly, Yp = r^x + Vyi). For 

simplicity we have neglected the source singularity. 

The effect of the axionic coupling is only seen when 
there are inhomogeneities in the material. Adding planar 
layers is identical to finding generalized reflection coeffi¬ 
cients except now we replace the source Eg with a vector 
containing the dyads. Thus the Green’s function for lay¬ 
ered TSB-TI’s is given by 


Gz^z'{r,r',u;) - j d^k, 

X m)!-' 


C(r,E) : F^^^,(z,z') 


h ^2 

fX'Zrip 


gikp-(rp-rp) 


with 


C(r,r') = 


m* (r) G m(r') n* (r) 0 m(r') 
m* (r) (g) n(r') n*(r) (g n(r') 


(87) 


( 88 ) 


and the : operator implying the element-wise Frobenius 
inner product. For z and z' in the same layer E^^^/(z, z') 
is given by Eq. ( [7^ and z') by Eq. (73). Eor 

z G n and z' G m in the different layers E^^^,(z, z') 


( f^ and F^^^,{z,z') by Eq. (|80|. 
reck, one can show that the result- 


is given by Eq. 

As a consistency c 
ing Green’s function reduces to that for a traditional 
magneto-dielectric material when the axion coupling van¬ 
ishes and that it satisfies the Schwarz reflection principle, 
which is required for the response to be causal (see Ap¬ 
pendix 1^ . 


VII. DIPOLE FIELDS CLOSE TO A TSB-TI 
SURFACE 

As an example of the use of the Green’s function we will 
compute the electric field pattern of a single frequency, 
dipole point source close to a TSB-TI surface at z = 0. 
We take the source to be in the upper layer, z' > 0. 


T(r,r',a;) z < 0, 

(89) 

each of which can be computed separately. Each part 
can be found by expanding the definition of the Green’s 
function given in Eq. (87). These expressions are given 
in Appendix [B} 


A. z-Orientated-Dipole 


First we will consider a single frequency, dipole point 
source, orientated in the z direction, placed close to a 
material surface. This source can be represented by a 
current density of the form 

J(r') = —iujd{uj)S(r')z, (90) 


where d{uj) is the dipole strength. The source is placed in 
the upper layer, which is taken to be the vacuum (e = 1 
and fi = 1), at 1.5 Ag above a surface, where Ag is the 
vacuum wavelength. The material parameters for the 
surface are e = 16 and ^ = 1, which are comparable to 
those of Bi2Se3 [36]. Substituting the current source into 
the expression for the electric field in Eq. (83) shows that 
the relevant components of the Green’s function are the 
G*^(r,r',w), where i = x,y,z depending on the desired 
field component at r. The expression for the Green’s 
function components can be simplified by converting to 
polar coordinates, after which the angular integral can be 
performed analytically. The resulting Hankel transform 
integral, however, must be computed numerically (the 
relevant integrals can be found in Appendix [C| . 

The field patterns for this configuration are shown in 
Eigurej^ Eiguresj^ (a), (b) and (c) show the, x, y and z 
components, respectively, for the real part of the electric 
field (equivalent to the time dependent field at t = 0) in 
the X — z plane for 0 = 0 - the case of a conventional 
magneto-dielectric. In this case the mixing coefficients 
vanish and hence one sees y-component of the field is 
zero. The field patterns for the x and z components are 
those that one would expect from a point dipole. Eig- 
ures|^(d), (e) and (f) show the, x, y and z components, 
respectively, for the real part of the electric field in the 
X — z plane for the case of a TSB-TI with, 0 = tt. In 
this case the mixing coefficients are non-zero. The axion 
coupling causes a rotation of the polarization of the field, 
generating a non-zero y-component at the interface. Eig- 
ures 111(g), (h) and (i) show the, x, y and z components. 















10 



FIG. 8: (Color online) The field pattern (real part) in the x — z plane for a single frequency, z orientated, point dipole close to 
a TSB-TI interface. Here, the upper layer is the vacuum (^ = 1, e = 1 and 0 = 0) and lower layer is a medium with e = 16 
and p = 1. (a), (b) and (c) are the x, y and 2 components, respectively, when the axion coupling in the medium is 0 = 0. (d), 
(e) and (f) are the x, y and 2 components, respectively, when the axion coupling in the medium is 0 = tt. (g), (h) and (i) are 
the X, y and 2 components, respectively, when the axion coupling in the medium is 0 = —tt. All distance are in terms of the 
vacuum wavelength and all field strength units are arbitrary. Note that, for clarity, the amplitude of the y component of the 
field has been scaled by x 1000 compared to the x and 2 components. 


respectively, for the real part of the electric field in the 
X — z plane for the case of a TSB-TI with, 0 = —tt. This 
case similar to that of a a TSB-TI with 0 = tt, except 
that the interface causes the field polarization to be ro¬ 
tated in the opposite direction. Thus, the y-component 
is the opposite of that in Figure[^(e). Note that the dis¬ 
continuity at the interface in Figures (c), (f) and (i) is 
expected since this is the longitudinal component of the 
electric field, which, unlike the transverse components, is 
not continuous at the boundary. In fact, as the difference 
in permittivity at the interface is 16, one would expect an 


order of magnitude jump in the longitudinal component, 
which is observed. 


B. x-Orientated-Dipole 

Next we will consider a dipole source orientated in the 
X direction, which is given by a current density of the 
form 

J(rO = - 


iu}d{uj)6{r')x, 


(91) 
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FIG. 9: (Color online) The field pattern (real part) in the x — z plane for a single frequency, x orientated, point dipole close to 
a TSB-TI interface. Here, the upper layer is the vacuum (^ = 1, e = 1 and 0 = 0) and lower layer is a medium with e = 16 
and iJ, = 1. (a), (b) and (c) are the x, y and 2 components, respectively, when the axion coupling in the medium is 0 = 0. (d), 
(e) and (f) are the x, y and 2 components, respectively, when the axion coupling in the medium is 0 = tt. (g), (h) and (i) are 
the X, y and 2 components, respectively, when the axion coupling in the medium is 0 = —tt. All distance are in terms of the 
vacuum wavelength and all field strength units are in arbitrary. Note that, for clarity, the amplitude of the y component of the 
field has been scaled by x 1000 compared to the x and 2 components. 


where, again, d(w) is the dipole strength. We consider a 
similar geometry as before with the source was placed in 
the upper layer (vacuum) at 1.5 Aq above a surface. The 
material parameters for the surface are, again, e = 16 
and ^ = 1. Substitution of the source current into the 
expression for the electric field in Eq. (83) shows that the 
relevant components of the Green’s function, in this case, 
are the G“(r,r',a;), where i G x,y,z depending on the 
desired field component at r. Again, by converting to po¬ 
lar coordinates and computing the angular integral one 
arrives at a Hankel transform integral that must be com¬ 


puted numerically. The relevant integrals can be found 
in Appendix |D| 

The field patterns for this configuration are shown in 
Figurej^ Figuresj^ (a), (b) and (c) show the, x, y and 2 
components, respectively, for the real (t = 0) part of the 
electric field in the x—z plane for a conventional magneto¬ 
dielectric (0 = 0). Here, as with the z-orientated dipole, 
the mixing coefficients vanish and hence one sees no y- 
component to the field. Figures (d), (e) and (f) show 
the, X, y and z components, respectively, for the real 
part of the electric field in the x — z plane for the case 
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of a TSB-TI with, 0 = tt. One, again, sees the gener¬ 
ation of a non-zero y-component owing to the effects of 
the interface. Figures (g), (h) and (i) show the, x, y 
and z components, respectively, for the real part of the 
electric field in the x — z plane for the case of a TSB-TI 
with, 0 = —TT. As before we see the inversion of the 
y-component compared to that in Figure ^(e). Again, 
the discontinuity in Figuresj^ (c), (f) and ^ is expected 
since this is the longitudinal component of the electric 
field. 


VIII. SUMMARY 


We have constructed the Green’s function of a layered 
TSB-TI and used it to study the field pattern of a single 
frequency point dipole close to the surface of a topolog¬ 
ical insulator. Reflection and transmission from a TSB- 
TI surface leads to mixing of TE and TM polarization 
components and hence a rotation in the overall polariza¬ 
tion of the incident light. This effect has the potential 
to be the basis for a number of novel optical and quan¬ 
tum optical effects, for whose study the Green’s function 
will be useful. Owing to the ubiquitous nature of the 
Green’s function in both classical and quantum electro¬ 
magnetism, it is hoped that the closed form expressions 
for this function will prove to be beneficial to a great 
many fields. 
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Setting uj —lo* gives 

Vx ^ VxG(r,r',-a;*) 

y{r,-uj*) 

+ i— — [V0(r, -uj*) X G(r, r', -w*)] 

C TT 

- ^^-Q-e{r, -a;*)G(r, r', -lo*) = (5(r - r'). (A3) 

As the permittivity, e(r,a;), and permeability, /r(r,a;), 
via causality arguments, also obey the Schwarz reflection 
principle and the axion coupling, 0(r, w), is real, we have 


Vx ^^^VxG(r,r',-a;*) 

-I- j— - [V0(r, w) X G(r,r',-a;*)] 

C TT 

- ^ £*(r,a;)G(r,r^-a;*) = S{r - r^). (A4) 


By comparing Eq. (A4) with the complex conjugate of 
Eq. 


_|), one can see that the Schwarz reflection princi¬ 
ple holds for axionic materials. 

For purely imaginary frequencies one has 


G*(r, r', i^) = G(r, r', - = G(r, r', i^), (A5) 


where ^ is real. Hence at imaginary frequencies the 
Green’s function is real. In addition one has. 


^- = \l^ (A6) 

where Hz is real. Substituting this into the expression 
for the Green’s function in Appendix converting to 
angular coordinates, performing the angular integration 
and noting that the resulting Hankel transform is real 
one can see that the components of the planar half space 
Green’s function obey the Schwarz reflection principle. 


Appendix A: Schwarz Reflection Principle 

The Schwarz reflection principle states that 

G*(r,r', a;) = G(r,r', (AI) 

One can show that the same principle holds for axionic 
materials. Given the Helmholtz equation in Eq. ([^, the 
definition of the Green’s function reads 

V X —-V X G(r,r',w) 

— i—— [V0(r,w) X G(r,r',a;)] 

C TT 

- ^e(r, w)G(r,r',w) = 5(r - r'). (A2) 


Appendix B: Greens Function for a Planar Half 
Space 

The half space Green’s function has a single interface 
at z = 0, hence d± = 0. Hence, the multiple reflection co¬ 
efficients = J. For a source point in the upper layer, 
the upper reflection and transmission coefficients vanish, 
,R_i_ = T-|. = 0, and the lower reflection and transmission 
coefficients becomes that for a single interface, R_ = R 
and T_ = T- Similarly, for a source point in the lower 
layer, the lower reflection and transmission coefficients 
vanish, R_ = r_ = 0, and the upper reflection coeffi¬ 
cient becomes that for a single interface, Rj^ = R and 
= T. Putting this together, expanding the dyads, 
and ignoring the extracted singularity, leads to expres¬ 
sions for the reflective, transmissive and free space parts 
of the Green’s function for a source and field point in an 
arbitrary layer. 
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The free space part of the Green’s function is propor¬ 
tional to J and is given by the well known expression 


Go(r,r',w) = /i(w) 


V G V -I- J 


oikR 


47ri?’ 


(Bl) 


R^^kp,k,,z,z') = 

/ /\ kxkz „ ky 

— sgn (z ) RtM,TM + -j^RTE,TM 


, (B8) 


which, on evaluating the derivatives, becomes 


G'o(r,r',w) = ^(w) 


oikR r 


47ri? 


. ikR — 1 \ , 

2z?2 


-k- 


3 — 3ikR — k'^R^ R 0 R 


fc2i?2 


i?2 


(B2) 


with R = r — r'. 

The reflective part of the Green’s function can be writ¬ 
ten as 

G(r,r')= [ ^e^^<’^’^-’^'^^R^Hky,kz,z,z'), (B3) 


(2^)2 


with 


R^^{kp,kz,z,z') = 


2kz 

y r> n 

^Rte,te - ^Rtm. 

Ivp n^piv 


TM 


k 

-l-sgn (z') {kxkyRTE,TM — kxkyRTM,TE) 

KpfZ 


(B4) 


Ryy{kp,kz,z,z') = 




(Id+b'l) 


2kz 


h.2 u2j^2 

j^UtE^TE — 7 2^2 ^TM,TM 

Kp /vp/C 

k 

—sgn {z') {kxkyRTE,TM — kxkyRrM.TE) 


h‘^h 

nipi\i 


(B5) 


i?"^(kp,fc„z,z') = 


2kz ' 


jfe^(|z|-|-|z'|) 


k h h h h‘^ 

•^x^y '^x'^y^z D 

-HtE.TE - 7 2 7 2 ^TM,TM 


U2 


h2h2 

rVpf^ 


k 

-l-sgn (z') (kyRTE,TM + klRTM,TE) 

fCp/v 


Zkz 


, (B6) 


k h k h 

^x^y T-) ^x^y^z D 

-KtE.TE - 7 2 7 2 


U2 

Up 


h2h2 

hjp Tv 


TM 


k 

— sgn (z') {klRxEXM + kyRTM,TE) 

fCp Kr 


, (B7) 


i?""(kp,fc„z,z') = ^e*'=^(l^l+l^'l) 

sgn (z') -^j^Rtm,tm + -^Rtm,te 


Zkz 

/ / \ ^y^z j-f kx j~. 

—Sgn [z ) —^Rtm.tm — ^Rte,tm 


/ f\ D 

) ^2 Rtm,tm — i^Rtmxe 


, (B9) 


(BIO) 


, (Bll) 


kp 

j^Rtm,tm 


(B12) 

where is the permeability in the layer. 

Similarly, the transmissive part of the Green’s function 
can be written as 

G(r,r')= J ^e^'^Eir.-<^T^^{kp,kz,z,z'), (B13) 

with 

r’^’^(kp, fcz, Z, z') = (‘^) ^i/c^|z|+zfc,/|z'| 

'^V rr I '^xEzEz' rji 

Pi^TE,TE + 1tM,TM 

fXip rvp fx rx 

1 f k > k 

-l-sgn (z') p- ikxky-^TxEXM + kxky-^TxMXE 

(B14) 


T^^(kj,, Z, z') = ^ikz\z\+ik,,\z'\ 

kf.kzkz' 

^^te,te + 1tm,tm 

IX p rxprxrx 

1 f k ' k 

—sgn(z') — [kxky-^TTE,TM + kxky-^TxMXE 

(B15) 
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T"^(kp,fc,,z,z') = 


2kz' 


ikz \z\-\-ik^i \z' I 


k* k k* k / 

rri ^X^y^Z^Z' rri 

~^1te,te + ^ 1tm,tm 

njp rVp A/ rv 

+sgn (z') — (^y-^TTE,TM — kl-^TTM,TE 


(B16) 


T'^^{\ip,kz,Z,z') = ^ik^\z\+ik^,\z'\ 

^ k* k* k k k / 

rri ^X^y^Z^Z' rji 

--j^lTE,TE + 1 tM,TM 

hjp njphjhj 

— Sgn(z') (^I-^Tte^TM — ky^TTM,TE 


(B17) 


/ f\ kxkz ^ ky ^ 

sgn (z ) Ttm,tm + -^Tte,tm 


, (B18) 


T^^{kp,kz,z,z') = y^Me^k,\z\+^k^,\z’\ 

j f\ kxkz' ^ ky 

sgn [Z ) Ttm,TM + -j^TTM,TE 


, (B19) 


Ty\\^p,kz,z,z') = ^^e^kM+^kA^'\ 

sgn [z') Ttm,tm — -^Tte,tm 


, (B20) 


simplified by converting to polar coordinate, after which 
the angular integral can be computed analytically. The 
resulting Hankel transforms, which must be evaluated 
numerically, are 

47r 7 

^2 

X Ji {kpRp) -^Rtm.tm, (Cl) 




1 

47r 


^2 

X Ji (kpRp) -j^RTE,TM, (C2) 


47r J 


k^ 

X Jo (kpRp) Rtm,tm, (C3) 


for the reflective part and 


X Ji (kpRp) ^ Ttm,tm, (C4) 

KK Kz^ 


Ty%r,r') = ^ J 


X Ji (kpRp) Tte,tm, (C5) 


T^y{\ 


■pjkz,: 


z') 


ifcjz|+»fc,,|z'| 

2kz> 


X 


Sgn (z') 


kykz' 

kk' 


Ttm,tm 


Ttm,te , 


(B21) 


r^^(kp,fc^,z,z') 


ikjz\+ik,,\z'\ 

2kz' 


kk' 


Ttm,tm 


(B22) 

where /i'(w) is the permeability in the source layer. 
Note that if the cross-reflection, Rte,tm and Rtm,te 
and cross-transmission, Ttm,te and Ttm,te, coefficients 
vanish one recovers the Green’s function for a conven¬ 
tional magneto-dielectric [52]. 


T^^(r,r') = £: y 

X Jo {kpRp) J’ Ttm,tm, (C6) 

KrC 

for the transmissive part, where Rp — Yp — v'p. 

Appendix D: Green’s Function Components for an 
x-Orientated Dipole 

The field from a x-orientated dipole close to a TSB-TI 
can, again, be found from the planar half space Green’s 
function in Appendix Converting to polar coordinate 
and evaluating the angular integral leads to 


Appendix C: Green’s Function Components for an 
z-Orientated Dipole 

The field from a z-orientated dipole close to a TSB-TI 
can be found from the planar half space Green’s func¬ 
tion in Appendix The relevant components can be 


r k 

X I [Jo (kpRp) + 2J2 (kpRp)] ^Rte,te 

k„k 


— [Jo (kpRp) — 2J 2 (kpRp)] -j^Rtm,tm ]■, (Dl) 
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X \ [</o (kpRp) — 2J2 (kpRp)] -^Rte,tm 
+ [>^0 (kpRp) + 2J2 {kpRp)] -^Rtm,te \ , (D 2 ) 


Ty^(r,r') = -^ [ dfcpe*'=^l^l+*'=xN'l 
Stt J 


X I [Jo (kpRp) — 2J2 (kpRp)] -j^TTE,TM 

- [Jo {kpRp) + 2 J 2 {kpRp)] ^Ttm,te\ , (D5) 


R^^(y,v') =-/ 

Att j 


^2 

X Ji (kpRp) j^Rtm,tm 


for the reflective part and 


T^“=(r,r') = — / 

Stt ./ 


(D3) 


X I [Jo (kpRp) + 2J2 (fcpi?p)] -j^Tte,te 

+ [Jo (kpRp) — 2 J 2 (fcpi?p)] , (D4) 


T^“=(r,r') =-^ J 




X Ji (kpRp) ^ , (D6) 


for the transmissive part. 
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